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The sliding Luttinger liquid approach is applied to study fractional topological insulators (FTIs). 

We show that FTI is the low energy fixed point of the theory for realistic spin-orbit and electron- 
electron interactions. We find that the topological phase pertains in the presence of interaction that 
breaks the spin invariance and its boundaries are even extended by those terms. Finally we show 
that one dimensional chiral anomaly in the Luttinger liquid leads to the emergence of topological 
Chern-Simons terms in the effective gauge theory of the FTI state. 


I. INTRODUCTION 

Topological insulators (TIs) are materials that exhibit 
a spectral gap in the bulk and at the same time have 
gapless excitations on their surface. The most famous 
example of a TI is the quantum Hall effect, observed in 
two dimensional conductors subjected to a perpendicu¬ 
lar magnetic field at particular filling fractions 1 ®. Re¬ 
cently, another type of TIs, symmetric under time rever¬ 
sal (TR) became a focus of experimental and theoretical 
research®). The experimental observation of these TR 
symmetric TIs has been reported in two dimensional 5 7 
and three dimensional material^' 

For non-interacting fermionic TIs and topological su¬ 
perconductors (TSs), the full classification based on 
their symmetry class and spatial dimension has been 
developed- 1315 . Among TIs, the simplest class of 
fermionic insulators invariant under time reversal T is 
the All (or symplectic) class, which is described by a Z 2 
topological invariant in 2+1 dimensions. A topologically 
non-trivial state within this class possess a helical edge 
statd^. 

For interacting systems the classification of TIs- 17 1 ^, as 
well as their microscopic description, is a subject of an 
active research. A toy model construction belonging to 
this class is a system of two non-interacting quantum Hall 
bars subjected to opposite magnetic fields, placed one on 
top of the other. The ground state is then composed 
of two species of electrons with opposite spin polariza¬ 
tions. As was pointed out in Ref [T^l in a properly tuned 
state an inclusion of electron-electron interactions within 
each layer leads to two copies of fractional quantum Hall 
states with opposite fillings. This induces the formation 
of a fractional TI (FTI) stat^22MH ; stable with respect to 
time reversal symmetric perturbation^ 1 ^. Although this 
model correctly captures the topological properties of the 
system, several simplifying assumptions were made: (a) 
The interactions between electrons in the different layers 
in the bulk was neglected, but taken into account between 
the gapless edge modes, (b) The spin-orbit interaction 
was assumed to be S z conserving and was chosen to cor¬ 
respond to a constant spin-dependent “magnetic held”. 

A convenient way of describing TIs and to account 


for interactions is the sliding Luttinger liquid (LL) 
approactPSHS Within this framework interactions be¬ 
tween parallel Luttinger liquids are engineered to de¬ 
scribe different nontrivial phases in two dimensions. De¬ 
spite its somewhat artificial appearance, the model cor¬ 
rectly captures the topological properties of different sys¬ 
tems and correctly reproduces the tenfold classification 
of noninteracting Hamiltonianj2Sl. 

Connection with previous works.- Within the sliding 
LL approach, the construction of the integer and frac¬ 
tional TIs was carried out in Ref. ;2B] and in Ref. E3 
In Ref. [26] the authors construct the different classes 
of topological systems based on antiunitary symmetries. 
This was done without stating any microscopic realiza¬ 
tion of the Hamiltonian that realizes the symmetries in 
question. In Ref. [27iKlinovaja and Tserkovnyak suggest, 
for the first time in this context of sliding LL approach, a 
Hamiltonian with spin orbit coupling that possesses TR 
symmetry. 

In that work the following assumptions were made: 
a) spin-orbit interaction was described by Rashba term, 
with the S z component only, b) The interaction between 
the electrons was assumed to be sufficiently strong to 
stabilize the desired fixed point, but was not analyzed 
explicitly. 

In our work we use sliding LL approach, explicitly 
taking into account generic electron-electron interaction 
terms consistent with TR symmetry in a model micro¬ 
scopic Hamiltonian. We also allow for a spin-orbit in¬ 
teraction that acts on all the components perpendicular 
to the wire direction (z and y ). We analyze the sta¬ 
bility of the FTI fixed point by deriving and explicitly 
solving the corresponding RG equations. We use this 
model to obtain an effective Chern-Simon theory as the 
low energy description and show that it correctly cap¬ 
tures the low energy physics of Abelian FTI states. Al¬ 
though the connection between Chern-Simons theory and 
the edge modes of topological systems is well knowrU^SI, 
our work is the first that links the sliding LL approach 
with a Chern-Simons low energy description. 

This paper is organized as follows: In the first section 
we review the wire construction of fractional (Abelian) 
quantum Hall states. We then extend the analysis to ac- 
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count for realistic Rashba and electron interactions. In 
the second section, we study the relevance of the multi¬ 
particle hopping operators that drive the system into the 
topological nontrivial state. We find that the interlayer 
interaction makes these operators more relevant, in com¬ 
parison with the toy model limit. This generic spin or¬ 
bit interaction is limited by the condition that it does 
not close the gap, which would lead to a transition to a 
different state, (e.g. through the appearance of a non¬ 
trivial spin texture). Finally, we discuss the emergence 
of the low energy description that follows from gauge 
invariance in the problem. Integrating out the massive 
modes we show that the low energy model is indeed given 
by abelian Chern-Simons theory, as expected for a FTI 
state. 


II. FTI FROM COUPLED QUANTUM WIRES 
A. Laughlin states from Luttinger Liquids 


After bosonization, the sliding LL Hamiltonian "Hsll = 
Ho + "Hfs takes the general (still quadratic) form 

H S ll = [ dx ( d x<f>j) TM jk(dx<f>k), (4) 

n jk J 

where (f>J = (<pj , 6j). The 2x2 forward scattering matrix 
is Mjk = Mij+Vjfc, with Yjk parameterizing the forward 
scattering interactions. At the filling fraction 

v = 2k°/b = 1/m, (5) 

momentum conservation allows for the construction of an 
infinite set of inter-wire many-particle tunneling opera¬ 
tors without fast oscillating terms (Friedel oscillations). 
Among these operators, the most relevant - in terms of 
renormalization group (RG) analysis - is of the form 

Oj = exp [ipj - ip j+ 1 + m(dj + 6j+i)\, ( 6 ) 


In this section we briefly review the coupled wire con¬ 
struction developed in Refs [33 and 132] for classifying 
topological states of two dimensional electrons. This ap¬ 
proach was recently applied for anomalous quantum Hall 
effeclP^, the Halperin states in the FQHEP^ an d for the 
construction of non-Abelian states in FThpm2H. We start 
with an array of parallel identical uncoupled wires sep¬ 
arated by a distance d. In each wire we place the same 
density n e = k°/ird of spinless fermions with single par¬ 
ticle dispersion E(k). These fermions are subject to a 
magnetic field perpendicular to the plane formed by the 
wires. We choose the gauge A = —Byx for the vec¬ 
tor potential. The magnetic field shifts the Fermi mo¬ 
mentum in each wire by the amount Skp j = bj, with 
b = \e\dB/h. The linearized low energy Hamiltonian, 
around the Fermi momentum kpj = yk 0 + bj for each 
chirality (77 = ( R,L ) = (+, —)) is 


Ho=v° f ^2 ( dxrji/l 3 {-id x - ( 1 ) 

jiV 

Under bosonizatior l 35 * 36 l (see A 1), this Hamiltonian be¬ 


comes 


= dx W j f + ( 2 ) 

i 

where dQj fit and tpj are the density and phase fields at 
wire j respectively. Two body density interactions are 
described by the Hamiltonian 

Hfs = ^Y, /( 3 ) 

jk,T)r]' 


which hops electrons between wires j and j + 1. In pres¬ 
ence of this tunneling term it is convenient to define the 
so called link fields^ 


2v? J+ ! = ipj + ipj +1 + m(0j - 0j+i), (7) 

2^+i = <pj - v j+ 1 + m{0j + 9 j+ i). (8) 


Using the commutation relations between the density and 
phase fields (A2) it’s easy to see that the links fields sat¬ 
isfy [d X <pe(x),ipt'(x')\ = [d x 9e(x),6t>(x')] = 0 and 


[d x Qi(x),ye{x'j] = inm5wS(x - a;')- (9) 

In the link fields basis, the Hamiltonian becomes 

H = H S ll + Y, [ dxgcos{28 t ). (10) 

e d 

When the cosine term is relevant under RG analysis, it 
opens a gap in the spectrum. This operator can always 
be made relevant by putting an appropriate choice of a 
forward scattering interaction. In this state, the system 
possess an excitation gap and quasiparticles characteriz¬ 
ing a Laughlin state at filling v = 1/m. 


B. FTI for odd m integers 

To construct the TI state we consider spin-orbit inter¬ 
action in each wires. The most general non interacting 
Hamiltonian on each wire j, quadratic in momentum in¬ 
cluding Rashba terms, confining potential and transla¬ 
tional invariance along the direction of the wire (Fig. [l]) 
is 
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FIG. 1. Wire arrangement considered for the construction 
of fractional Hall states and FTIs with time reversal. In the 
former case a constant magnetic field Bz is assumed. For the 
FTI construction, the presence of Rashba interaction and a 
special external confining potential is assumed (see text). 


Hj = +a so (p x B) ■ VV(y,z) + V(y,z) 

2m e 

= 2m - + (A* 0 * + ^j a v)P* + Vj, (11) 

where m e is the electron’s effective mass, aso is the 
strength of the Rashba coupling and Vj = V (jd, 0) is 
a spatially dependent confining potential at the position 
of each wire (here we assume that the plane formed by 
the wires is located at z = 0). The parameters Aj’ z are 
simply 



FIG. 2. Dispersion relations Ej(k) for the different wires. Due 
to the Rashba term, the energies of different eigenstates are 
displaced to opposite sides around zero momentum. The ever 
increasing Rashba terms induce a different momentum shift 
for each wire. The arrows in the top of each parabola indicate 
the corresponding spin projection. Around the Fermi level, 
tunneling processes are allowed (see main text) and they are 
represented by red (green) arrows for the + (—) spin states. 


Note that each Kramers pair is defined up to a phase. 
This signals the explicit the break of SU{ 2) invariance 
(due to the spin-orbit coupling) to a 17(1) x 17(1) sym¬ 
metry. 

To reach the topological phase (see also the discussion 
after eq. (401), we tune the confining potentials such that 


W dV 

X J= aso ^ 


y—jd 1 
z = 0 


dV 

X z = -a SOir 

J oy 


y=jd' 
z =0 


( 12 ) 



and V 2 



The simplest potential V that leads to the topologi¬ 
cal nontrivial phase corresponds to a parabolic confining 
potential V(y, z ) = V\y 2 /2 + v 2 yz. This potential gener¬ 
ates a space dependent Rashba coupling with strengths 
Xj = —(cxsovid)j and Aj = ( aso‘V 2 d)j (compare with 
Ref. IHSland Eq. (3) in Ref. f27|l . 

The dispersion relation Ej{k) for the Hamiltonian (11) 
is (with h = 1) 


Ej(k) = ^ n + ^ ± ky/( A~) 2 + (A")2, (13) 


The eigenstates of (11) are 


i>j,+ ( x ) = e 

v jr (.r) = e ik * 


_ ^ikx ( ^ S1H Oij 

cos a a 


COS Oij 

—i sin OLj 


(14) 

(15) 


with tan 2ay = — yt • Due to the time reversal symmetry 


of the Hamiltonian (111, Tipj s is also an eigenstate with 


the same energy (Kramers partners). 


The energy dispersion becomes in this case E ± (k) = 
(7’±jA) 2 /2m e . In the case of a single quantum wire, this 
energy dispersion can be obtained also by a combination 
of spin orbit and a Zeeman field. In such scenario it has 
been showrP 7 ! that the system possess fractional quan¬ 
tized conductance and can host Majorana bound states 
by proximity-coupling with a superconductor. 

With the choice (161, the Fermi momentum in the wire 
j becomes 


k F, s U) =Vk° ~ sAj, (17) 

where A j = j A and y = ( R,L ) = (+,—) denotes chiral¬ 
ity. In the previous formulas the values of y = R/L are 
understood as ±1. The values of s = (+, —) corresponds 
to the different spinors. For convenience we will refer to 
electrons with s = 1 (—1) as those that belong to the 
upper (lower) layer. This choice of potential results in 
a particularly simple dependence of the Fermi momen¬ 
tum on the wire index. This will enable us to construct 
relevant multi-particle tunneling operators that conserve 
momentum and are free of oscillations, analogously to 
how it was done in the Section Hi Al 
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Using_^j a as a basis, we proceed to linearize the Hamil¬ 
tonian Jnl around the Fermi energy. The linearized 
Hamiltonian is 


H 0 = v° f J2 j dx{rf R s (-id x - k Fs {j))tpf s 

■F, s {j))' t Pj,s}- ( 18 ) 

General density-density interactions between wires are 
given by 

HfS=1 ^ H / dx PjA X ) V vv',ss'Pj', S '( X )- ( 19 ) 
jk,rjrj' ,ss' 

where pj s (x) = ipjs^js is the density of electrons in the 
wire j , belonging to the layer s, moving in the direction 
rj = ( R,L ). Using the stan dard bosonization rules de¬ 
scribed in the appendix (A4) the bosonized Hamiltonian 
reads 


is local in terms of the new fields (f)^ = 


2(p 

2 9 


o+i 


0+i 


= 


= fs 


v: 


j+i 


<7; 


J+ 1 


s\m\(0l-9l +1 ), (25) 

s\m\(0l+9i +1 ). (26) 


Following Ref. j^Uwe name these bosonic degrees of free¬ 
dom link fields (note that now they have an extra layer 
index s). From the definition of the link fields and the 
commutation relations of the bosonic fields 9 and ip, de¬ 
fined in (A2), we find the link fields’ commutation rela¬ 
tions 


[<pe,s{x),9e',s>{x')\ = wrs|m|sgn(:r - x')5 aa >Su>- (27) 

The quasiparticle charge density and current on the 
link i and layer s are defined accordingly 


n = H 0 + n FS = f dxidJYMijidxfi), ( 20 ) 

i,j 

where 4>j is the vector of bosonic fields at the position j, 
c/S = $+, i(P_, 9 3 _) T that are associated with fermions 

from the upper (tp R + ,4’f,+ ) and the lower (ip R _,ipk_) 
layers. 

At this stage, the 17(1) x 17(1) symmetry, which ap¬ 
peared after the explicit break of 517(2) by the spin or¬ 
bit term, becomes manifest under bosonization as the 
symmetry under a constant shift of the phase fields 
ip{x) —► <p(x) + /?. 

The most general quadratic Hamiltonian with nearest 
wires interactions has the following structure 


•0 


sdxO^s 

jQ,e,s 

= PQ/,s = 

7r|m| 

•1 


sd T 0i, s 

Jq,£,s 

= JQAs = 

7r|m| 


and (28) 

(29) 


Note that the link fields defined above transform under 
time reversal as 


Tpe, s T =-ipe,- a + 


(1 — S)7T 


T9, s T~ l = -0,-s 


(30) 

The 17(1) x 17(1) symmetry of the original fermions 
becomes a symmetry under constant shifts in the links 
phase field ipi, s . This is no longer the case when the sym¬ 
metry is gauged. This leads to interesting consequences 
for the low energy theory, described in the next sections. 


Mij = MgSij + j. (21) Using (30) one can show that the most general form 


In order to represent the Hamiltonian in terms of lo¬ 
cal link variables, cf. Eq. ([7]), we choose the forward 
scattering matrix of the following form 


of the forward scattering matrix M, consistent with time 
reversal symmetry, is given by 


( 22 ) 

(23) 


M 0 = UHlU + MV, 
Ml = MU, 

with U and V given by 


U = 


For this special forward scattering matrix (21) the 
Hamiltonian 

R = ^2 f dx(d x <j>t) T M(d x (f> e ), (24) 

* J 


M = 


^ OL(p(p OL<pQ CXipip \ 

OiipQ OLQQ Oi^Q OLqq 


Olu 


x <p8 


a v e 


(31) 


\ a v e a ee a <p0 a se ) 


/i 

m 

0 

0 V 


( i 

—m 

0 

0 \ 


m 

0 

0 1 

!■ -1 

1 -i 

m 

0 

0 

0 

0 

1 

—m 

, 0 

0 

1 

m 1 

Vo 

0 

1 

—m ) 


^ 0 

0 

-1 

—m ) 


being parametrized by six independent variables (a vv> , 
a^g, OL V (p, a^g, a g g, age). In terms of the link fields the 
Euclidean action is given by the sum of over all the link 
fields 


S 0 = ^2 J dxdTi(d T (f) e ) T K(d x (f) i ) - V.. 


(32) 


Here K is the matrix 


I< = 


27t|to| 


a x 0 
0 —cr x 


(33) 
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After integrating out the bosonic phase fields (pc, S: we 
are left with the action for the fields 9c, s describing den¬ 
sity fluctuation on the links. This action can be cast in 
a familiar form by using the charge and spin fields 


0q+ = 


t,+ + Qe,- 

V2 ' 


0 /.- = 


i,+ 


V2 


(34) 


In terms of Qc,± the action becomes (summation over 
a = +, — and links l is implied) 


Sn = — [ dxdr 
2tt ./ 


Ufy , ,, \2 1 


JT&W + K u 


(d T ®e, a ) 2 


+ic((d x Qc, + )(d T Q e ,.) + (d r 0*, + )(d x 0*,_))J(35) 
with Luttinger liquid parameters K± given by 


1 a ee ± a gS (a v e ± oi s ) 2 

—— (36) 


q/2 — q2 

tpip l P t 4 - 5 


The velocities u± in terms of the interaction parameters 
become 


‘Z'KTTlX / T Q(pip) 


aoo zb a 6 Q - 


(q y e ± a vS ) 2 

Ot^ptp A OC^p(p 


while the parameter c (which explicitly breaks parity) 
reads 




27r|m|(c 


) 

5 ) ' 


(37) 


The solution of this equation exists for all values of j 
only if Aj is linear in j. This restricts us to the confining 
potentials given by Eq. (16). 

Setting 


A j = j A 


we rewrite Eq.(40) as 


ak° _ Y(pP( S p,a + S p,a) 

~K~ ~ 


E 


qR _ qL ^ 

p\°p,a °p,a) 


(41) 


(42) 


R s L are in- 


where we have used that a = (+, —). As s 

.0 

tegers, equation (42) has solutions only when is a 
rational number. The simplest FTI phase corresponds to 


ak° 

~A 


a 

2m 


(43) 


This is similar to the condition ([5]) in the construction 
of Laughlin states. In that case, if the filling fraction 
matches 1/m, the tunneling operators ([6]) conserve the 
momentum and can be relevant. In analogy with the sin¬ 
gle layer scenario, eq. (42) corresponds to the condition 
of an effective filling fraction i/ a = 2ak° /A in the layer a. 


Any operator 

ters {Sp a ,Sp a }, satisfying Eq.(42) describes a legitimate 


with the set of parame- 


multi-particle hopping. The most relevant operator in 
the RG sense is given by 


The action (351 is self dual under p -<->■ a and resembles 
the action that appears in the context of spin ladders 1 ^. 


1. Tunneling operators 


,R 

am — 1 

Sl,a = 

am + 1 

1 ,a 

2 

2 

,R 

am + 1 

,L 

am — 1 

0 ,a 

2 

b 0,a 

2 


s R / a L = 0 for 0,1 


(44) 

(45) 

(46) 


So far we did not allow for electron tunneling between 
different wires. To account for such processes we consider 
the multi-particle hopping operators 


/-ql s p,<j> s p,a} 

3,a 


n tyf+pA x )y*' a tyj+pA x )y*' a - ( 38 ) 


which are integers when m = A/2 k° is an odd integer. 

{ s R ,s L } 

Upon bosonization, O' a p ’“’ p ’ in this case becomes 

O t , s = cos(20/, a ), (47) 


Here we use the convention if -1 = ipf The possible 
choices of integers s R and s L are restricted by charge 
conservation E p ( s p,a + s p,a ) = 0, and momentum con¬ 
servation 

J2 S P,a k F,aU +p) + Sp,a k F,aU + P) = 0, (39) 

V 

replacing the value of the Fermi momentum we have 

k ° J2( S P,a ~ S p,a ) - A 3+p( S p,a + S^J = 0. (40) 

p p 


where 9c is defined in Eq. (25). The corresponding action 
is given by 


S = S 0 + gJ 2 

i 


dxdr (cos (29c, + ) + cos(2 Be-)) ■ (48) 


As we see, in the bosonic notation, the inclusion of tun¬ 
neling operator leads to the sine-Gordon type action. In 
the case when the cosine term is a relevant perturbation, 
a gap opens in excitation spectrum of the sliding LL in 
the bulk. To study this question we realize an (RG) anal¬ 
ysis of this operator. 
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III. RG ANALYSIS, THE BULK GAP OF FTI 
AND EDGE MODES 

A. Relevance of tunneling operators 


The RG equations can be derived for this problem in 
the standard wav^SI, and one finds 


d_ 

Jl 


K, 


d 

Jl 


‘Up-h'P 




d 

Jl 

d 

dl V u„K, 


U(j 

Ua 


1 


Ag 2 f(x) 


A g 2 f(xf 


(49) 

(50) 

(51) 


Here x = u p /u a is the ratio of sound velocities in the 
charge and spin sector. The functions f(x) and h(x) are 


f{x) = \Jh(x, c) + Jh(x,c) 2 - x 2 , (52) 

h(x, c) = (1 + x 2 + |c| 2 a;)/2. (53) 

The scaling dimension of the tunneling operator is 
given by 


A[cos(26b, s )] = 


K 0 + K a 


2 \ /1 + 


\c\ 2 x 

(1+ai) 2 


(54) 


There is no renormalization of the parameter c up to 
order g 2 . 

The set of equations ( 49pl ) can be linearized around 
the fixed point (A 0 , <jr) = (2,0). The fixed point 
A(K p , K a , x) = 2 defines a region (in parameter space) 


Kp + K a 


= 1 + 


(l + ;r) 2 ’ 


(55) 


for the LL parameters and velocities. Expanding around 
the point p° = (K pl Kjx°) belonging to the surface 
defined in Eq. ((551) , we define A = 2 + A with A <C 1 


and y = g/u* (see also appendix (A 4)). The set of RG 
equations (valid in the vicinity of the fixed point), is given 
by 

dy 
dl 
dX 

Jl 

where the terms quadratic in A were neglected. The pa¬ 
rameter C = C(K p , K°, x°) is an involved function of the 
point p° = (K° Kjx°) which belongs to the surface 


Ay, 

(56) 

Cy 2 , 

(57) 


(551. 


For C < 0, the cosine term in Eq. (47) is never rel¬ 


evant and the system goes into the weak coupling fixed 
point. In a generic situation C > 0 and the tunneling 
operators become relevant for A < 0. This corresponds 
to a Berezinsky-Kosterlitz-Thouless (BKT) transition at 
A < 2 (see also (A23)). 



FIG. 3. Renormalization flows for Eqs. (56 and 571. The dif¬ 
ferent lines indicate the different initial strengths of g/u p . The 
thick dashed line represents the separatrix of the BKT transi¬ 
tion. When the initial g/u p is above g* /u p , the system flows 
to the strong coupling regime (fractional topological phase), 
while for initial values below g* /u p , the system flows to weak 
coupling regime (sliding LL). This plot is made for c = 0.1. 


As shown in Fig. [3j this implies that the system flows 
to strong coupling regimes if 


I<p + K a 


< 2W1 + 


(1 + a;) 2 


(58) 


For c = 0, the theory becomes massive when the LL 
parameters K p + K a < 4. For generic interactions (pre¬ 
serving time reversal symmetry), c / 0 and the region 
where the Luttinger parameters flow to strong coupling 
is enhanced as shown in Fig [4] In this sense, generic time 
reversal symmetric interactions help driving the system 
into the topological phase. 

From now on, we assume that the system is in the 
massive phase. Our RG eqs. show that the cosine term 


in Eq. (47) is relevant and a gap in the spectrum of 


sliding LL opens (this was just assumed in Ref. [27]). For 
periodic boundary conditions, the gap opens everywhere, 
as the field = (p™ — s\m\9^ in the wire N pairs up 
with the field (J s = ip° s + s\m\6° s in the first wire to form 
the link fields, all of which develop a gap after the cosine 
term becomes RG relevant (see eq. (251). 


B. Edge modes 

Experimentally however, a two dimensional finite size 
sample is modeled by open boundary conditions. In this 
case note that the first and last fields (with different 
“layer” components) 


= ¥>° + s\m\e°, (59) 

Vn = <Ps ~ s\m\6 f, (60) 
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(a) 



(b) 



FIG. 4. Various separatrix lines with various interaction 
strengths c. (a) For a generic electron interaction (finite c) 
the critical point A shifts to the right with increasing c. (b) 
The slope of the separatrix lines decreases as a function of c, 
that means that the region where the system flows to strong 
coupling (i.e the topological phase) is enlarged. 


chiral states appearing in the edge of a disk, an array of 
concentric wires can be used, as discussed in the appendix 

El 


Few comments are in order here. First, we note that 
our analysis is valid only for special filling fractions, for¬ 
mally described by Eq. (42). Moreover it is limited by the 
special type of spin-orbit interaction in Eq. (41), and the 
confining potential in Eq. (16). As we are about to see, 
in this case the system has topological order. Although 
we do not address the stability of this point, we expect 
that static disorder stabilizes the topological phase, mak¬ 
ing it robust against small deviations on the filling factor 
and generic spin-orbit interactions. 

Our analysis is limited to the Laughlin states with an 
effective filling factor v = ±1 jm , to being an odd integer. 
A heedful reader could ask, what does prohibit us from 
repeating this procedure for even values of to. While 
formal steps of construction remain valid, the result is 
not the true ground state. As we know from the context 
of fractional quantum Hall effeclP^ the true ground state 
for these filling is highly sensitive to the fine interaction 
details, and thus requires special consideration, that is 
beyond the current work. 

We now show how our approach can be connected with 
a more conventional Chern-Simons theory. 


IV. LOW ENERGY DESCRIPTION AND 
CHERN-SIMONS TERM 

Until now, we have considered the dynamics of 
fermionic excitations in the system. We will now dis¬ 
cuss the electromagnetic response of the system and the 
effective low energy description of the FTI state. As ex¬ 
pected for the TI case, the resulting theory indeed has a 
topologically non-trivial Chern-Simons term. 

To find out the electromagnetic response, one needs to 
compute the effective action I as follows, 


had no cosine term and therefore remain ungapped. 
These fields at the boundary correspond to the gap¬ 
less edge modes of a TI. The same conclusion has been 
reached in Ref. EH At each edge we have two counter 
propagating edge modes, labeled by s = (+, —). The dy¬ 
namics for these edge modes is described by the Hamil¬ 
tonian 


exp(-/[A]) = J 


(62) 


Here S is the fermionic action in the presence of an exter¬ 
nal gauge field A. The coupling of the gauge field to the 
fermionic action in its bosonic form in Eq. (48) should 


be gauge invariant. We proceed to investigate the gauge 
invariance of S below. 


Hedge = ^ m j VF y f dx{d x (/) s 0 ) 2 + {d x (j) s N ) 2 . (61) 

These helical edge modes are the hallmark of the topo¬ 
logically nontrivial phase. We will see in the next sec¬ 
tion that the phase that we encountered here is described 
in the long wavelength limit by a Chern-Simons theory, 
which makes the topological properties more evident. 

Different boundary conditions can be implemented us¬ 
ing different arrangements of wires. In order to reproduce 


A. Gauge invariance 

Let’s recall the gauge transformation laws for the gauge 
field A „ and original fermions ’>pj(x. t). Due to charge 
conservation, the phase of electron operators can be rede¬ 
fined i/j'(x,t) = e la ip(x,T ), without changing the action. 
This (7(1) freedom can be gauged by making the phase a 
(initially a constant) a function of space-time. Of course, 
this change alone does not preserve the fermionic action. 












To make the action invariant again it is necessary to in¬ 
clude a gauge field with transformation laws that can¬ 
cel the extra contribution from the space-time dependent 
phase. We have then the following gauge transformations 

tyj(x,r)] G = (63) 

[ A nd( x , r )] G = A hA x > t ) - d^aj(x,r). 

It is easy to check that these transformations leave in¬ 
variant the free electron action coupled to an external 
gauge field 

S hee [ip 1 A\ = j cPx'tpiid^ - A^^ip, (64) 


i.e. S llee [ip G , A G ] = S iree [ip,A], (here ip = ipA° and 7 0 ’ 1 
are the Dirac matrices defined in (A3)). 

Bosonizing the electron (see ( |A ![ )), we can read off the 
gauge transformation rules for the phase ip and density 
field 6 


[<pI{ x ,t)] g =(p J s (x,T) + a j (x,T), (65) 

[°IA,t)] G = 01(x,t). 


The link fields in Eq. (25), that appear as the natural 
basis after the inclusion of the tunneling operators (47), 
change non trivially under a gauge transformations 


the tunneling operator for the s = + (upper layer) fields 
in terms of the original fermions is given by 


oi,+ = ^u) ' f^,,.) ^ ' , 

(67) 

and similarly for the lower layer (with to replaced by 
—to). Under the gauge transformations (63) the tunnel¬ 
ing term On changes as 


it 


[O j+i , s ] G = O j+hs e^- a ^\ (68) 

In order to maintain gauge invariance, the operator 
Oe tS should be modified, introducing a Wilson line Wjj+i 


W jd+ 1 = exp 



such that the tunneling operator Cb iS becomes 


(69) 


O 


j+b s 


o j+hs w jJ+1 . 


(70) 


It is clear that this combination is invariant under a 
gauge transformation, as the phase acquired by O is ex¬ 
actly canceled by the transformation of W 


Wj+\A x ' t )] G = <^+i,s( a; ’ r ) 

[0 j+ i, 5 (a;,T)] G = 6 j+ i iS (x,T) 


Otj{x,T) + Otj + l{x, t) 
2 

a j{x, T ) - 07 + 1(35 i) 
2 


Note in particular that s transforms as a gauge field 
with ctj — ay+i the discretized version of the derivative 
along y. 

The gauge field parallel to the links, along with the 
temporal component (scalar potential) couples to the 
action via the current and density operators. After 
bosonization the part of the action that couples the gauge 
field to the quasiparticles is explicitly 


d2x ie,s A K = 


'2 


e.s 


d x 




( 66 ) 

where is the quasiparticle density (n = 0 ) and cur¬ 
rent (/i = 1 ), defined in (28); is the two dimen¬ 
sional Levi Civita antisymmetric tensor, with e 01 = 1 
and summation over repeated indices is assumed. The 
gauge field along the wires is A\\ = (Aq(t,x),Ai(t,x)) 
(/i = ( 0 , 1 ) = (r, x) in | 66 l» 

The last term in Eq. ( 66 ) can be integrated by parts, 
and ensures that the action is gauge invariant with re¬ 
spect to the temporal and the x-component of the gauge 
field. The gauge field A x = A 2 which is perpendicular to 
the wires induces an Aharonov-Bohm phase in the inter¬ 
wire tunneling operator Oe defined in Eq. (47). Indeed, 



The bosonized action (including the modified tunnel¬ 
ing term (70)) becomes 



This action is gauge invariant (here A\ = A y j and we 
have used that J A y dy = A y d , where d is the distance 
between wires which is assumed smaller than any other 
length scale). 

When the cosine term becomes relevant, the field 0e 
is pinned to the minimum of the potential, i.e 9e = 
—Afd/2. The action acquires a term proportional to the 
chiral anomaly in each effective layer (+,—), and this 
anomaly dominates the dynamics of the system after the 
massive modes (8 fields) have been integrated out 39 l In 
this case the effective action for the upper layer becomes 
(using Ay = d) 
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I+[A] =2d^Aid,A {W + .... (73) 


where ... denote terms containing higher derivatives, 
that do not contribute to the small momentum behav¬ 
ior. Here we recognize the discretized version of (see Eq. 
(A 5) for more details) 


I+[A] = 


1 


47t| 7 


d?xe^ p Apd„A p , 


(74) 


that is the Chern-Simons theory for the fractional quan¬ 
tum Hall effect with filling fraction v = 1/m. The effec¬ 
tive action I + accounts for the low energy dynamics of 
the system, describing the electromagnetic response with 
respect to the gauge field A p in the upper layer (similar 
construction can be done for the lower layer separately). 
Taking the functional derivative respect to the gauge held 
A p , one expectedly reproduces the Hall conductance 


H + [A] _ y. _ 

5 A 


fi 27t|to 


e^ p d v A p , 


(75) 


where ,7/ denotes the electron current in the upper layer. 
For the lower layer we have similarly 


I~[A\ = —^ [ d?xe^Apd v A p . (76) 

47r|m| J 


To satisfy (identically) the conservation equation dp = 
0, it is customarto define the bosonic 6+ held by J p = 
A. e ^pQ v \)+ (note that these auxiliary bosonic helds b p 


are essentially the full gauge invariant version of e IJI 'd u 9(). 
The held theory in terms of 6+ that reproduces (75) is 
therP^l 


£ b [6±, A] = - b-d v b~) 


47T 


--e^Ap(dX + d„b~). 


2n 


(77) 


where we have included the contribution from both lay¬ 
ers. Note that the theory (77) can be cast in the form of 
the BF theorjEHHl, by defining the helds a p = bp + b~ 
and bp = b+ — bp. This double Abelian Chern-Simons 

theory U( l) m x U( l) m accounts for the ground state de¬ 
generacy m 9 of the system if placed on a genus g sur¬ 
face. This theory also de scribes t he exchange statistics 
of quasiparticle excitation d 28 ! 43 ! 44 ! . 


V. SUMMARY AND OUTLOOK 

Starting from an array of parallel LLs with Rashba 
type spin-orbit interaction, we have constructed a FTI. 
The magnitude of spin-orbit interaction depends on the 
position of the wire and mimics the magnetic held in the 
quantum Hall effect. The interaction between electrons 
is restricted to nearest and the next nearest wire, but 
includes all terms consistent with T symmetry. We have 
derived and analyzed the RG equations for the multi- 
particle tunneling operators and showed that there exist 
a window in the parameter space where these operators 
are relevant, and FTI is a stable hxed point. Due to 
the large number of possible interaction constants the 
precise position of this window can be determined only 
numerically. We have analyzed the stability of the FTI 
hxed point in relation to a different interaction constants. 
Remarkably, the stability is enhanced by inclusion the 
repulsive interaction between electron with opposite spin. 

To establish the connection with other methods, we 
have derived an effective low energy theory for our model. 
The resulting hxed point of our construction is captured 
by a double Chern-Simons theory with correct counting 
of topological ground state degeneracy and exchange of 
quasi-particle excitations, and agrees with one suggested 
in Ref ll9l within a more restrictive model. 

Before ending the paper, we list some of the questions 
that remain open: 

a) The stability of the FTI against perturbations that 
break TR symmetry. Though a strong perturbation of 
this kind certainly destroys the FTI state, the topolog¬ 
ical protection may hold against ’’weak” perturbations. 
For the magnetic impurities located solely on the edge, 
the BKT transition between TI and trivial insulator was 
founcP®. The case of magnetic impurities added in the 
bulk remains to be studied. 

b) Within the sliding LL approach the tunneling oper¬ 
ators are chosen for the states exactly at fractional (or in¬ 
teger) fillings. In the presence of TR preserving disorder 
one expects the state to form a plateau. The emergence 
of plateaus yet remains to be shown within sliding LL 
approach. While for non-interacting case the formation 
of plateaus and the transition between different filling 
is described within cx-model approach^, the current for¬ 
malism provides a natural framework in the presence of 
interactions. 
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Appendix A 


1. Bosonization 


Under bosonization, a collection of the fermionic operators labeled by index j, of chirality = (R, L ) = (+, —) and 
spin s = (t,l) = (+, —) around the Fermi point kp S - becom^HM! 


V’i.sOM) = 


TTJ 

- O ’ J 


72 


nx c 


(Al) 


The bosonic operators 0 3 (x,t ) and (p 3 s (x,t) follow equal time commutation relations 


[^ 0 , t), 9 J S , {x', t)} = iTTSgn(x - x')5 ss '5 j3 ’, 


(A2) 
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B 



and t// s is a Klein factor that ensure anticommutation of fermions. In this notation, chiral fields cf> J R s (right mover) 
and <j? L s (left mover) are 


4>r, s = 4 + 

The bosonization of electron densities become 


and K,s = e i- 


Pis = (4441 = o ( d x d l + d M) and Pu = (444* = - d x ipl). 


(A3) 


(A4) 


2. Set of concentric loops and tunneling operators 


In this appendix we show how the wire construction can be done for a disc geometry. We choose a set of N w 
concentric wires (see Figf]h with linearly increasing radius (i.e. radius of circle j being Rj = jRo). We classify the 
states by integer n, the component of angular momentum component in the direction perpendicular to the disc. The 
geometry of the problem makes it convenient to use the he symmetric gauge A(r) = ^f-0. Solving the Schrodinger 
equation for spinless electrons one finds energy levels for the wire j 




h 2 / n . eBR^ \ 2 

2 M e R 2 \ j~ J 2 ch J 


(A5) 


It is clear that this set of wires also coincide with the positions of guiding centers of the wave functions for the 
same gauge, Landau level and angular momentum. The sliding LL construction can be perceived either as specially 
fabricated setup, or as a complete basis of highly anisotropic states. 

For a given chemical potential // = 2 m r 2 1 Fermi (angular) momentum in the wire j is ripj = ±jn^+j 2 &(Ro), 

with Q(R) = fh/^) c the magnetic flux in units of the flux quantum. Note that ufj is well defined just for $(I?o) € Z. 
We are interested in the tunneling terms of the form 


Oj,j+i = (4 J+ i) O’rj+i) " ¥ ‘ "' 2 * 


in terms of bosonic degrees of freedom ip, 0 the tunneling operator (A6) becomes 


Ojj+i = 


_ i(n° F m-$)(2j+l)x 


exp (<pj - (fi j+ 1 + m(0j + ©j+i)). 


(A6) 


(A7) 


At filling fractions v = n° F /' I> = 1/m, the phase of this operator cancels. For appropriate interactions, the operator 
Ojj + 1 becomes relevant under RG flow. In this case, all the modes but the last one are gapped, and the system 
becomes topologically equivalent to a quantum Hall annulus. This results in the single chiral state at the circumference 
of the disc. The similar construction for TIs takes into account the degeneracy between Kramer’s partners, and is 
done using a set of concentric wires hosting spinfull electrons. 
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3. Time Reversal Invariance 

Time reversal is implemented by an antiunitary operator T with T 2 = — 1 for spin 1/2 particles. We use the 
representation 


T = lT X IJylC, (A 8 ) 

where t x acts in the chirality index as t x ^ )S = and cr y Pauli matrix acts on the spin index by <J y ipri,s = 

(cx^) Sjs'^Prjjs' ■ represents complex conjugation. 

For (l+l)-dimensional topological insulators, we can study the classification of time reversal systems by examining 
the representative Dirac Hamiltonian for spin 1/2 particles 

H = J dxi[) s (—i^dx + m)ip s . (A9) 


with = ip'^o- The two component spinor i/j s contains the chiral fields for each spin s. The matrices 7 1 and 7 0 
satisfy { 7 ^, 7 ^} = g^ u - We use the chiral representation 


7 


0 



(A10) 


Introducing a regularization scheme (studying the system on a iV-site lattice for example) the Hamiltonian (in first 
quantization) becomes a 2N x 2 N block diagonal matrix, where each N x N block acts nontrivially in the chirality 
subspace. The time reversal invariance of ( |A9| ) is reflected in 

THT~ X = H -> ( ia y )H T (-ia y ) = H, (All) 


where we have used the hermiticity of the Hamiltonian. Note that the last equality in (All) is the definition of the 
symplectic Lie algebra sp(2iV). Hamiltonians satisfying (All) are said to belong to the All or symplectic class. The 
evolution operator exp (iHt) is then an element of the symmetric space (coset) U{2N) / Sp(2N). 

If we bosonize the fermionic Hamiltonian (A9), we obtain 


1 + /» 

H = 2 7 T ^ 7 dx(d x 9 s ) 2 + ( d x ip s ) 2 + — cos(20 s ) 


(A12) 


Using (A 8 ) and (Al) we find that the bosonic fields transform under time reversal as (here s = —s) 


T9 S T 1 = 9 S , and Tip s T 1 =-ip s + 


1 — s 


(A13) 


so the Hamiltonian (A12) written in terms of bosonic fields is invariant under TR (as it should be!). Note that 
although (A12) is a non-quadratic Hamiltonian in the bosonic fields, still we can associate to it a symmetric space, 
inherited from the fermionic description. 


a. Special points in parameter space 


It is worth noticing that some regions of the parameter space in the interacting Hamiltonian (20) correspond to a 


non interacting system in terms of link fields. In particular the region parametrized by = a v g = a e § = 0 and 

(aggacpcp — a g) 


= I<~ 2 = K~ z = K 


a 2 w m 2 


-2 


—2 


A'r: 2 m 2 (agga tfitp - a 2 g) = u p = u a = u. 


(A14) 


corresponds to two independent Luttinger liquids with the same Luttinger parameter K. At K = 1, the link fields 
can be mapped to non interacting fermions. 
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4. Renormalization of Tunneling operators 


The RG equations (491 can be linearized around the fixed point g = 0 and A = 2. This last point defines a surface 
in parameter space 


K p + K c 


= 1 + 


(1 + x) 2 ’ 


(A15) 


for the Luttinger liquid parameters and velocities (x = u p /u a ). Around the point p° = (AT 0 , K °, x°) lying in the 
surface (A15), we expand A 

A(p + <Sp) ~ 2 + + ( SK P + 6K a ) — + a . 0 )[(i + a;0)2 + |c|2a:0 ) 5x (A16) 

where <fp = (i 5K p , 5K ai 5x ) and |<5p| <C 1. The RG equations (49) linearized around the point p and g = 0 become 


dg 

dl 


K° P + RT p 


(< SK P + SK a ) - 


|c | 2 (1 — x °) 


(1 + £°)((1 + X 0 ) 2 + |c| 2 x°) 


5x 


|(^ p ) = -[/°(1 + (/°) 2 )( A ' 0 ) 2 ] (^) 


dl 


(6K a ) = - 


f 


,((x 0 ) 2 + (/ 0 ) 2 ) 


0 \ 2 




!( fa >= - {ff) - K ° Mx ° )2 - (I)' 


(A17) 

(A18) 

(A19) 

(A20) 


where /° = f{x°), dehned in (52). Defining the variable A = K ol K o (SK p + SK a ) - j| e | 2 a; o ) Sx we can 

study the renormalization of A and y = g/u° p . Combining the last three equations we get the BKT type of RG 
equations 


^ =-A y and ^ = -C(K°,K°,u°,ul)y 2 


controlled by the sign of the function C. Defining the always positive function 


A = 


2/°((l + (/°) 2 )( K) 2 + ((z 0 ) 2 + (f°) 2 )(K 3 ) 2 /x°) 

K° p + K 


(A21) 


(A22) 


C is explicitly 


/ 4c 2 (l-a: 0 )i 0 

' (l-f 2 )K° + (f 2 -(x°) 2 )K°Jx 0 ' 

^ + (l + cr°) 3 (R° + A'°) 

(1 + / 2 )(AT 0 ) 2 + (*0 + f 2 / x°){K®) 2 ~ 


(A23) 


where the values (AT 0 , A'°, a; 0 ) are related by the equation A(AT° K°,x°) = 2. Note that for x -C 1, x 1 and x = 1 
C is positive. 


5. Discretized Chern-Simons and continuous limit 


The actions (73) and (74) differ in the bulk by the term / d 3 x (A x d y A T — A T d y A x ). This term appears 
disguised in the wire construction. It is related to gapless edge modes in the x direction in the construction based on 
Luttinger liquids, for a finite size system. These gapless edge modes in the x direction are inevitably connected to 
the seemingly missing term 47r j m j f d 3 x (A x d y A T — A T d y A x ), as is shown below. 

Let us recall how the bulk Chern-Simons term is related to the edge modes as a consequence of gauge invariancd^ 47 -! 
The action for the hydrodynamic field 6 + is given by (771 
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CS +[b} = M| d 3 xe^b+d„b+ 


(A24) 


where we take the effective layer + for simplicity. Under a gauge transformation [ 6 +] G = 6 + — the gauge 
transformed action C<S' + [ 6 ] becomes 

[ C5 +[ 6 ]]G = M ^ = M ^ ^^(6+ _ a M a)a„6+ 

i.e [C75 + [6 ]] g = CS ,+ [&] - ^ [ d 3 xe pup d p ad v b+ (A25) 

47T 


where we have used that e pvp d u d p a = 0. The last term in (A25) can be integrated by parts in the manifold 
E = Ex {[0, -C/a:] x [0 ,L y ]}, leading to 


[C\S ,+ [ 6 ]] g -C , S ,+ [&] = 


47T 


dr / dya (d y b x - id T b+) 


' —oo J 0 


C=I/rr 


V ' ItC —0 


/ dr dxa (d x b x — id T b x ) 
l — oo J 0 


\V~Ly 

\y-0 


(A26) 

assuming that the gauge transformation a vanishes in the infinite past and infinite future. It is evident that the 
Chern-Simons action alone i s not invariant under a gauge transformation in a manifold with boundary. One way to 
solve this problem (making (A26) vanish identically) is to impose that the gauge field 6 + beco mes a pure gauge on 
the boundary, i.e. b p \dT, = o^tp. this resolution makes the fields at the boundary dynamical as (A24) becomes 


CS + [b] = j dr I dx j dy {b+ {d x b+ - d y b+ ) + 6 + (<9 y &+ - id T b\\ ) + 6 + (id T 6 + - d x b+) } , (A27) 


’ —oo J 0 J 0 

=Mru‘'4\ 

’-oo Jo Jo 


+ 


2 tt 

M 

47T 


Km - M) +aiA'UUUAfo) 


/ 


nOO 


/ dr dx(b+b +) 
l—oo J 0 


I y — Ly 

ly =0 


/ dr dy{b+b+) 
' —oo J 0 


— L x 

y v r j ix=o 


(A28) 

(A29) 


Using the condition of pure gauge at the boundary in (A29) we have 

b+{id T b+) -bt\id T b+) 


CS+[b] = 5 i\ b r^ b v- d V b t) + 


+ ^ { I dr I dx {d x (t) + ){id T (j) + )\ v Lv - I dr I dy {d y (/) + )(id T (j) + )\ x 


— L x 


47T 


' —OO j 0 


y=o 


' —oo J 0 


I 21 = 0 


(A30) 

(A31) 


note that the boundary term (d x (f> + )(id T (l) + )\ y= Q V appears integrating by parts b x d y b T — b T d y b x . In the partition 


function Z = f Vbe cs ^ the hydrodynamic field b+ is integrated over. The path integral over b T imposes the 
constrai nt (d x b+ — d y b x ) = 0 which makes the fields b x ,b y pure gauge in the whole manifold. This implies that the 
r.h.s of ( A30| vanishes. The Chern-Simons action depends just on the fields in the boundary of the manifold 


CS+[&] = ^ 

L J 47T 


/ dr dx (d x (j) + ){id T <j) + ) 
' — oo J 0 


+ \ | y—Ly 

y =o 


dr dy {d y (j) + ){id T (j) + )\ x 


' —oo J 0 


x—0 


(A32) 


Finally, to connect with the wire construction, we take L x — > oo, obtaining 

I I /* OO nOO 

C'5 ,+ [t>] = —— / dr dx [{d x 4>t){id T <t>i) - {d x <j>t){id T <t>t)] , (A33) 

which describes the chiral edge modes in the first and last wires (for the + layer). As we discussed, these edge modes 
appear as boundary terms after integrating f s (b+d y b+ — b+d y b+) by parts. Coupling the hydrodynamic field 6 + 
to an external gauge field A p and integrating over generates the term f s (A x d y A T — A T d y A x ). 

Including the action from both layers, and using L y = dN , with d the distance between the wires, we have 
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CS+[b]~CS-[b] = ^ J2 


dr / dx [(d x cj) a N ){id T (j) s N ) - (, d x (j) s 0 ){id T (j ) s 0 )]. 


(A34) 


s=+,— 


' — oo «/ —oo 


This action together with the Hamiltonian describing the gapless edges (61) amounts to the full description of the 
edge dynamics 


= Fl y 

47T ' ^ I . 

' s __l J — oo J — oo 


dr dx[(d x <t> s N )(sid T <j> s N +v F d x cl) s N ) + (d x <t> s Q )(-sid T cl)o + v F d x <l)o)]. (A35) 



